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ABSTRACT 


A detailed examination of the literature pertaining to subccooled 
critical heet flux was conducted, Subcooling, system pressure, ress 
flow rate, system geometry and flow stability were determined to be 
the important parameters in determining critical heat flux, The effect 
of each of these is discussed in detail. 


Various schemes for correlating experimental critical heat flux 
were subjected to critical examination, Ferm, content and comarison 
with experimental data were the bases for evaluation, The correlation 
of Griffith was found to compare most closely with the experirental 
data over a broad range of variables. A linearized cerrelating equation 
was developec fer the MIT low pressure (30 to 90 psia) data, including 
that determined as part of this study. 


Critical heat flux was determined experimentally for pressyres 
from 30 to 90 psia, mass flow ~ates of 3.69 x 10° to 11.76 x 10 
lbm/ft¢-hr and diameters of 0.047 to 0,2)2 inches ever the subcooled 
region. The experirental data obtained is in general agreement with 
that of other investigations. This low pressure data, when combined 
with the results of previous investigations, provides a fairly complete 
picture of subcooled critical heat flux for tubes up to 0.25 inch 
diameter, A graphical interpolation scheme is presented by which it 

is possible tc obtain predictions of critical heat flux hich are 
generally superior to that of the matheratical correlations. 


Thesis Supervisor: Arthur E. Bergles 
Assistant Professor of Mechanical "neineering 


$ 





TABLE OF COMPENTS 


TITIZ PAGE 
ABSTRACT 
TABLE OF CONTENTS 
LIST OF PIGURES : 
I. INPRODUCTION 
II, PROCEDURE 
TIT. RESULTS 
Factors Influencing Critical Heat Flux 
Correlating Squations | 
Results of Experimental Study 
IV. CONCLUSIONS AND RECOMMENTATIONS 
TMiscussion of Experimental Results 
Discussion and Evaluation of Correlaticns 
Ve APPENDIX 
A. Small Mameter Critical Heat Flux Map, 
Recommended Frocedure for Use 
Be Figures 
C. Deta 
D. Apparatus 
E, Yxperinentel Ffrocedure 
F, Sample oxperimental Data Sheet 
CG. TDevelormment of the MIT Small Diameter Correlation 
HK, Nomenclature 


I, eferences 


Nm 


I wv & w 


15 


33 


39 
Li 
53 
60 
67 
70 
71 
72 
7h 


aa 


may 





Q 





r 


ure ‘Yo, 


Ww NO 


LO 


LUST aT FTCGR S 


Scheratic Diagram of Test Avparatus 

Tetails of Test Section Assevbly 

Tependence of Critical Heat Flux on Subecoling 
forsD = ,.079 in., p = 30 psia 

Dependence of Critical Heat Flux on Subecoling 
for DP = ,.079 ine, p ® 90 psia 

Nependence of Critical Heat Flux on Tube Diameter 
Tependence of Critical Veat Flux on Subcooling 
for D = ,079 and .2h2 in., p 2 90 psia 

Dependence of Critical Meat Flux on Heated Leneth 
Comparison of Correlations and Exverimental Data 
for D> ,2h2 ine, p = 90 psia 

Compariscen of Correlations and Experimental Nate 
for 1D = ,180 in., p 2 2000 psia 

Dependence of Critical Neat Flux cn Subcooling fcr 
PD = .079 ing, p = 17, 347 and 735 psia 

Dependence cf Critical Tieat Flux on Subeccling for 
D® .079 ins, p = 1102, 1:70 and 180 psta 
Dependence of Critical Heat Flux on Subeccling for 


™ = .079 ines Dp ® 2205, 2570 and 290 psia 


Pace No, 
2 
he 


4,3 


Lb 
bs 


16 
47 


48 


Lg 


50 


51 


52 











» 
a 9-40 a a” a 


uy = on) hate SS & . an 


ar! Paras) 120 BHA) GO) Ge eat 






= Fg 28s. 2 a 

OP SE ee re, 6 Te! 
a ee ae 

oe SE Wig at Sw see 
a . DAG bee says 
Rae Cd, sawn 


ee = 
ies = ake re roat 7 
SERA ete 


a i} 
on Gites bac. “wo 


a7 





























































7 


a 
. 





ec ®& 







=— sas 4 
: 


oe 


L. INTRWELCTION 

In the continuing search for devices of hifher enerry flux Uensity, 
mcre and more enphasis is being pleced upon boiling as a means of 
increasing heat transfer rates with vroderate temporature iillerences, 
Current developments in high energy electromegnets, high verformance 
electron tubes and rcecket motors require very hich heat fluxes in 
cooling applications, Although heat fluxes well in excess of One 
Million Btu per square fcot per hour can oe obtained, heat fluxes 
camot be increased without limit. For a given set of flew ccnditions, 
there is a limiting Critical Heat Flux (CHF), at which the transition 
from nucleate to film boiling cecurs, At this ccint, the wall ternereture 
increases dramatically, usually resulting in catastrophic thermo-mechanical 
failure of the heat transfer surface. One also finds references in the 
literature to this as "burnout", and "the beiling crisis". 

Often, the most desirable mcde of belling fcr heat transfer 
applications is that oi surface, or subccooled boiling, where the average 
or bulk fluid temperature remains below that cf saturation. In this 
case, vapor formed at the surface condenses when it ccmes in contact 
with the ccoler main stream. The Pressurized Water Reactcr uses water 
at high pressure (2000 psi) as ccolant and ~cderator. The high pressure 
allows sufficient coclant terperatures for further heat transfer and 
power generation, while maintaining bulk temperatures below saturation, 
The maintenance of subcocling is important in reducing void fraetion 
and retaining continuity of the moderator, The electromagnets and 
electron tubes are generally cocled by low pressure systems where the 
expense and cemlications of high svstem pressures can be avoided. 
Utilization of subcocled boiling eliminates the complexities of two= 


phase flow, as well as the mechanical complication cf condensing 
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equipment; however, the CHF problem is et.il p esent in succocled 
boiling. 

Over the pect several years, mimerous stucies have been concered 
with this CHP vrobler, A great deal of experimental data is available, 
under various conditions cf flow. Investivseters have usually selected 
certain paramters and investirated the effect cf varvine these within 
certain limits, The majority o* the previous experimentation hac 
been conducted with unifcrm heat flux and round tubes. This investigation 
has been limited to this region, 

The first phase cf this investiration censisted cf a search of the 
literature to deterrive the flow ronditions which had been investirated. 
Techniques, generally in the form of correlating eeuations, heve been 
propcesed fer predicting CHF, by almest evervene who has ccnducted 
investigation in this field. One phase of this study has been concerned 
with a critical exarination cf the techniques and their icrm and 
agreement with available experimental data. 

Another phase cf this investigation ecnsisted of completion of 
Bome previous experimental work. Experiments were conducted at 
selected tube diameters, at low pressures and various flew rates, 

Date thus obtained, when combined with that in the literature, provides 
a picture cf CHF, fcr tubes cf small diameter and flow rates of interest, 
ever the ccmplete range cf prescures, from 30 osia tc just below critical 


pressure, 











<0 - 2 ; . =P ne ws ee £ 
frie WES 2 two part studys <¢ eWmmey, cl tee ciara (ec eecuure 


ane correlating editions anc exirimented Sctereinteien «fF “ri tceal 
Heat Flux (CHF) under certain os ii tions. 

First, citations of ali of the current literature concerned with 
Subecoled Criti@éal Heaw Flux were obleined. From these, the micrity 
of the literacure vas reviewed, en? experimental data an‘ oro ceed 
correlating equations were reccrie? fer further investircation. / 
large part of the work in this field has been ececmplished in Fussie, 
consequently certain information is available only in the Pussian 
Journals, untranslated. [I+ was agually pesvible to extract data and 
prepcesed ccrrelating equations without complete translation of the 
particular paper. 

A listing of propcesed correlating equaticns was prepared. These 
were compared with actual experimental data from varicus sources, 
including the experimental phase cf this study, in an attempt to 
deterrine the validity of the ecrrelastion schemes. 

Experiments were condurted to determine CHF for predeterminec 
values of the impcrtant variables. The dependence cf CEP uron 
subcocling was determined for three fixed values cf mass flow rate, 
with constant system pressure, tube length and diameter. This data 
completes the low pressure end of the stucyv cf Crmatski an? Kichigin 
(h, 5 and 6). Then, using cne of the same flow rates anc pressures, 
and helding exit subcocling constant, dependence cf CHF uson tube diameter, 
and upon tube length-diameter ratio was determined, Next, the variation 
of CHF with changing entrance forms, with the cther variables held 


constant, was investigated, Finally, tube diameter was changed te obtain 


compariscn data, anc to determine the maximum CHF available at the 





larger diameter with the exnoerimental opparatus, 


All of the experimental data relating te Subcooled Critical Mat 
t MIT was collectec, 


Flux that has been cbtaine? within recent yo 
converted to similar fcrmat and vlaced on machine data carts in 
anticipation of future studies. 

A detailed «description of the experimental apnaratus and prccedures 


emplcyec is contained in Appendix Tf and &, 
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FACTORS TW LUO" TG CRITICAY. MEAT FLUX 

Although there has been a , reat deal of research and isvestipetion 
in this fieldy particulary in the lor’ deeatn, there exists todar no 
clearly defined, universal, mathervatical descripticn of the precess 
which results iffa limiting CRP, Feol Sciling is itself a complex 
problem, the adiiticnal comlitiomcf subrcolec, forced fice tenting 
to cloud the picture even ore. Attemmts at dimensional analysis, 
am] treatment of the prcblem from basic egquaticns have, in tbe min, 
preduced little information, usually yielding leng, complicated 
functicnsl relationships end invelving simplifying assurpticns net 
necesserily borne out by experiment. 

It is generelly accepted, however, that. subcocling, system 
pressure, mass flew rate, flow cecmetry ant stability are the major 
factors in determining CHr. Thess variables evidence various decrees 
ef interdependence, and in discussing the effects of each, this 
interdependence rust be berne in mind, 

Effect cf Subccoling 

Subcooling is essentially a measure of the anount cf pre-heating 
required to cbtain boiling. This is generally measured as the 
difference between outlet bulk tenperature, or mere precisely enthalpy, 
and that at saturated ccnditicns, This is reasoratle fer unifern 
heating, in that failure normally cecurs quite near the damstream 
end cf the heated section, From ancther viewpcint, subeccline is a 
measure of the ccoling effect cf the main streerm. Thus CHF increases 
with subeccling. There is a tendency in the literatvre to treat CHY 


as linear with subcooling. This is not the case, but it is in ceneral, 


not a grossly inaccurate approximation, In the region of low values of 





subcooling the data show 3 


CRY is ~meliec. 
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Study, it. is st®ted that a variaticn of system prescule from ©) psie 

te WOOO psia ham little significant efiect on OFF: this in tubes ranging 
fror 0.3 inch to 0,09 inch in clameter. Hcwéver, tits is net ocrne 

out by thc clher expe imental data (9). A pressure increase trom 


600 pSia to 1000 msia results in a reduction of CHF ci as much as 25%, 


Further increases in pressure abcve 1900 psia result in further 


& 


* 


reductions in CHT, At lower precsures there is disagreement over the 
effect of presvurc, McAdams (32) and Gunther (33) repcrt nc eifect up 
to 163 psia, while Mirshak, et al. (31) shew a pressure dependence over 
the entire pressure range ccnsidered, Berges (11), reports a 15% 
increase in CH}, for an increase in vressure frer 20 to 85 psia. 
Chirkin and Iukin (2) demonstrate thet increasing pressure, in the 

Low pressure regicn, increases CHF. Their vressure effect is ruch 
stronger than that ci other investigators, however, it is now felt that 
most of their data suffere? fror instability prcblems which could have 
cicudee the real picture. 

The experimental data coes sugeest, however, that ior fixed geometry, 
tlew rete end subcocling under conditions cf stahle flow, there is some 
pressure fcr which CH attains a maximum value. The pressure at which 
this meximum value cf CHF occurs is net clearly defined, being dependent 
upon mass flow rate, and possibly the other varialles as well. It is safe 
to say theurh, that under conditions of reasonacle flow rete, this maximum 


eccurs scmeuhere between 100 and 600 psia. Unfortunately data in this 
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Vatiocac® of “i th Jianwetéer, Morceneic’: ae Tahiead (9), Pe jane 
tAis. eS kg 6 f. See Gl Meroe atm ea Ne wows fs Stl ve 
that the varliat@en of Cle Sith @xit subrccling is damier na LinAr, 
e8pecialiy at dbo vaye hig’ rresgeureas, FF is found, “hat @? the@ee high 


= 7m - ~ » -~ =. + a ¢ i 27 3 , a ~ 
pressures FTF becomes lest deeencent soem subcociing as smbecoliLag is 
& TI < ax 4. 4 x. ~ “a Bn pee woo” 7a t my 
inereased, This sugpests that at pressures noer cri ical, end with vo rv 


high values of sudcocling Abo ereazer then 200 “PJ/1b), CIT beceres 


Somewhat i-denerdent cr 


suoecoling. ‘ther hifh presvure tata shows 
this same general trend (1, 5 ar@ 4), This would, of crurse, require 


° 


further inves*sication tc be throurhly substantiate, 
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Generally, in the subcocle: region, ivirreaging mass slow ret 
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results in an increased (iY. This ic easily visusliszed since the 
increase? flew sweeps the wubble away from tie heat transfer surface 
mere rapid’y, replecing then with relatively cci¢e waver. Ecwever, the 


rate of ineresse in CIF cue tc increased Tlow rave is créss 
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At high pressures, the effect ef inereased mass Fler rate beccmes less 
prenounrec, 
E’fect cf Geometry 

Tire majority ef the research ond study ir thas fieid hes ben 
i‘watede tc rcun! tubes end arnular fluid paseages. The grees ei ects 
of the other variebles rewever, ere much the sare fer any cress section (8), 
In general, ter piven flow foil oreesure ane subeceling, CFi is greatest 


for e¢ cireuler tube are successively lover as Lhe geometry becomes rcre 


\ 


complex. ‘Swirl flow" is one emcepticn te this. The ratic oi heated 





perimeter to fotal wetted ne~imeter hes been cugpected bs ipg@ertant in 
the case of "swirl flew" (12). 

Iinvestisators disagree whetrer there is ¢ ram: where “lameler has 
no effect, Teroshchuck art Lantsman in a recent study (9), conduc ied 
tests fer foir tube diaweters over mos cf the high pressure rane, 
(735 psia to 2560 nsia), with exit ccnditions from high subeccling 
(Ah. 4.200 BTU/1b,,) to bulk boiling (X54 2002 de This ssudy showec 
CHF decreasing consistently with increasing channel diameter, iImrlicr 
investigations show decreasing CHF with decreasing diameter. sobth 
Chirkiin ear Tukin (2)) and Kafernhaus anc Becharov (26) report this. 
Again, this is now felt tc have been in a large part cue to flow 
instabilities in their systems, | TuFont (8), Mirshak (3) and 
Zenkevich (1;) note no diameter effect. Bergles (11) nctes a strong 
inverse dependence of CHF cn ciameter for small diameter tubes, but 
he shows the Zenkevich Prediction (1h) as a possible asymptote for 
Gameters freater than 0,3 inches. 

With flow rates high enough to eliminate the effects of natural 
convection within the passage, channel orientation has little effect 
on CHF, 

The effect of heated length has not been completely determined, 
In general, the length, or length-diameter ratio, is not important if 
the length is long, or L/D is large. The popular conception is that 
CHF 4s determined by purely local phenomena, thus as long as the flow 
is fully developed and entrance effects play no part, the length is 
not significant. But the question of how long is long enough has not 
been satisfactorily resclved, partly because of the intcracticn of the 
variables, in general, the Russian studies consider an L/D ratio of 


10 te be large encugh te eliminate the length effect, This conclusicn 
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has not been completly verified. Rereles (11), tn 9 Lei nressure study, 
shcwed an 1/7 dependence existing for L/7 on the order of 10. Doresheruch 
and Fried (18), nresent data fer tubes of 0,1? ineh Ciemeter at 1590) nsia 
that shewe nc differerce jn CHF for L/7 of 1b, 50 ard 1%, ever 2 fairly 
bracd rance cf exit ecnditions. ™Font (8) sugrests considers.’ on of 
the termera‘ ure grecient town “he heeted tube. The temrerylum: i erence 
is of ccurse dictated by the first Law of Thermodynamics, sc «ho. this 
pradient raally reflects is e lenethediametoer conside~stion. 

Tn general, CH! decreases for inereesine L/D, but fer hich oressures, 
the dependence cf CFF upen L/D becomes less sipnificant for values of 
L/P? greater than 30. 
Effect of Flow Stability 

This refers to vulsatiens in flew end vressure, resultine from 
cempressibility effects in th> heated sertions, If there is an expander 
cr accumulator upsteam cf the heated seeticn, there is a tistinct 
possibility of preducing a flew oscillation which will rericdically 
reduce the flow neer the tube cutlet end orrnaps cause a premature 
failure of the haot transfer surface. Bersles (11) noints cut that there 
are two trnes of instsbilttyvy: flow oscillation due to ccmmressiblility 
and flew excursion due tc system charactaristies. Flow conditions under 
which these pulsations de effert CHF are ecnsidered unstable. The result 
ef unstable flcw is a drestic reduction in CPP. with unstable flcw 
resulting ir heat fluxes as low as 20% of the stable velue, Stevie flee 
usually can be insured for subrooles boiling by installins some sort of 
throttling device (orifice, valwe, ete.) in the system just up=streanm of 
the heated seetion, Doreshchuck and Lantsmen used a cressure drop of 
1:70 psi in their high pressure studies (9, 19). Perrles (12) reported a 


pressure drop of appreximatel; 185 psi in his lew pressure work. 
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£. 05 12 ewi 08, Vhesstanilitvy evrch lee iq mesticular’ ecate in 
systers cottuainitig hest excrangers with maraliel iwoes from cowmcn 
inlet enc exit heacers, where a substaytiol emo «7 copmressl cle 
volume wey exist. In eny avert, bre ine all tian of © -westtfe ~oluciny 
feviee, ard trhe*resultant loss of downstream svetom pre sum, 36 more 


vr, 


than. justified by tae inerease in Om, 


Several other factors have been investicatec tc determine now chey 


Surface rcughnes’ , for reasonably emooth (commercial qelity) 
channels has little or no eiiect (7, 13). 

Tissolved pas in the flow strear, in szell ancunts, has a nergligibly 
erall efiert on CH (47), 

For beat flux unevenly distributed arcund the hested vccrimeter, 
peak heat flux was found te increase. For heat ilux unevenly distributed 
along the length of the heate? channel, quoiitatively similar results 
have been cbhtsined (7). Tests heave been conductes under ronditions where 
the heet flux was ferred into a cosine distribution along the tube and 


=_—— 


CFF was shovr to be reduced by scme 30% (11), 
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Numercus schetes have been vronosed as vorreédategy Cra mee, oe 
by which C¥F can suoposedly be predicted, The earliest c¥ trem ire 


of @ vewer function fcr: 


sé 


(c/b oy mAe EC TH at 


where B, " and F are the important variables, cr cethinaticns of them, 


(iD 


and c, @ and ¢ are experimentally ceter-ines exncnents. 


correleticns cf this tyre ten’ te neclet b4 in evrerticn botwveen 


t-4) 


variables; in cther eases the expcrents therstives orc cunctions co 
the controlling pavieweters. <Ancther bread categery inciudes the 


ne deta with 


i. 


Usuperprsiticn" cerrelations, those ce phbining veol dei} 


+ 


forced ccnvection heat transfer infermation, There cer-eletions are 
eften Jons anc difficult $c use. 
Several attes™pts have been made at anxlvytical treatment, csenkevich 
(11) and Pyabev and Lerzine (36) with dimersionel anslysis and Griffith 
23, 17) with a thecretical eralvsis., Cther equveticrc vere developed 
by “prese frree" date analysis, "sin> computiy -acrines, “aPent 


favers 2 form: 


where A, B and D are the influencing variehles, (pressvre, subeccling, 
Giameter, velerity,) and Cy Cys Coy ete, are constants, 

mm gemeral, the correlations are Jimited tc specific ranges of the 
variables, departing racically frem the experimentel vata in the excluded 
arees, in some eoses, cote finds ccrrelotions dervolcned Troe a pati cof 


mw e 
5 


= 
th 


an ‘exmeerimental study, weine certeivr date an’) oxchuling, “her. <C% 
techrique is usually accompanies bv a statement to tre efiert thet so 


many percent cf the experimental results agrees withir plus cr einus so much, 
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Tn cases such as this, the potentiel vror is chl’-en to examine the 
filterin: vrorecure <8 closely as he dces the correlstion itself, 

As a result cf this state of affairs, one phase cf this stucy 
was concerned with exemination of correlating equatiens. the listing 
of correlations is not all inclusive. Several ccrrelaticns were exartined 
and discarded, generally due tc their saikce similarity tc others. 

Of the correlation schemes exarined, the following are those 
which agreed reascnably well with the experimental data, or illustrete 
some point. 

Details of Farticular Equations 
Bell: (ref: 21) 


(a/Mor = 0.427 x 10° x 2.59 + 1.61 


“ (2,35) 1.77 “(s) 06 


Range of values: (a/A)ey = Btu/hreft°s p = 2009 psia3; G* .2x 10° - 


06 


5.0 x 10° Im/hr-ft2; Aga = 0 = 130 Btu/lbm 

This correlation is compared with experimental data (figure 9) 
for conditions very nearly identical, from three different sources, 
in obtaining this data from the various references, a certain variation 
was unavoidable; two people do not conduct tests at precisely the same 
conditions. However, the variation was kept to less than the expected 
experimental error, The data from reference 15 was obtained in 195h, 
prior to the general acceptance of the need for stability, and may be 
of questionable valus, 

It can be seen that the correlation follows the general trend of the 
data, with variations of 15% to 30%. 
Jens & Lottes: (ref: 17) 

This correlation was based on data taken at UCLA and Purdue over a 


ranse of pressures from 500 = 3000 psia. The values of K and n were 
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determined on the basis of the data such to make all deviations positive. 
(Q/A)on © K (c/10°)" +( Dtgat) a 
UCLA Data correlated t+ ~0%, +233, 0 = .26 ine 
i; ie gx. 
500 g—027 16 
1000 0626 = 4275 
2000 hiS = 50 
Purdue Data correlated to -0%, +605, D = .143 in. 
ce ae all. 
3.900 0915-9 275 
2000 655 50 
3000 630i 725 
The correlation was plotted for data (figure 9) ef appraximately 
D © .100 in,, G = 261 x 10° lom/nr-ft2, at p = 2000 psia using the 
constants and exponents of the Purdue data. 
Ryabov and Serzina: (ref: 36) 


1) (@/Moy = (a/A)np + 403 WeOAT IG! eg Toot 






Range of values: D= 2 mm = 079 ingy AT. & 16 = 590 deg I’, 

L/D = 20, p = 10 = 200 atmospheres = 147 = 2910 psia, G = 3.69 x 10° - 
22.4; x 10° Ibm/hr-rt? 

2) (a/Aep ® (G/M) gn * 257 Veo ATea¢eo? ky T,4°7 0515 





Range of values: D7 mm (.158 ine), Toat = 16 - 430 deg F, 

I/D = 10 = 220, V = 3 = 46 ft/sec, p = 76 = 210 atm © 1030 = 3090 psia 

Both of the above correlating equatims were developed from 
dimensionless equations. In each case (a/A) is the same and is obtained 


i? 


on 


wee @ 





from a plot of pool boiling CHEF derived as a function of pressure. hie 
correleting equatims were derived from a dimensionless equatio: 


O79 in.) 


1) Kyq = 3710 Kye® KypeO? Premed (9 = 20m 


iY 


or 2) Kaq * 2370 Ky 2 © K peor prvelt (D 5 hmm = .15¢ in.) 
where Ky, ® (Aq/A)om hep py J £6 


Ks Tgat (P21 “eve 





K ¢ = Alsat hee Py J So 


gv Tsat P2z 7? 
Pea 


It is assumed that (q/A).p = (a/A)., + (Aa/A)oy 

The 2 mm (.079 in.) correlation wes derived from the data ef Ornatski 
and Kichigin (1, 5, & 6) by the usual technique of letting the data 
determine the coefficients and expoents. Generally the date falls within 
- 30% of the correlating line, The correlatia was also computed for the 
2 mm date taken during the experimental part of this study, In this case 
the correlatio did not follow the basic trends of the data and was in 
error by as much as 190% at some subcoolings, 

The correlation for D = lh mm (.158 in.) was computed and plotted for 
Ge» 2.1 10° lbm/hr-ft*, p = 2000 psia and compared with the data for 
Deroshehuk and Lantsran's (9) ) mm and 8 mm tubes. The correlation fit 
the 8 mm (.316 in.) data very closely , but was 15% - 30% low on the 
h mm data. It was not plotted on figure 9 to avoid confusion as it 


fell generally om top of the correlation of Zenkevich and Subbotin. 
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Zenkevich: (ref: lh) 


ee Sa 





irate mph “ae 


(q/Aor * hee (233 & on eo) x [2 5 + 18h Site | x 107? 


Range of values: p * 1; - 165 psia, Ahgay = 20 = 262 Ltu/ibn, 

G = .22 x 10° ~ 9.96 x 10° Ibm/nr-ft2, D = 16 - 50 in. 

This iced a rs is linear with subcooling as commonly noted for 
Russtén analyses. Specific instructions as to the temperature to 
evaluate viscosity at are lacking. Using the fluid saturetion temp- 
erature yields an equation linear in subcooling, while using fluid bulk 
temperature prodiices a slight curvature in the resulting plot. This 
correlation is errr in figure 8 for both conditions and it can be seen 
that the equation dees not really reflect the experimental results. 
Zenkevich and Subbotin: (ref: 1) 

(a/A)er = 397 C2 AT gat *23 (ve/veg)te® — Btu/nrett? 

Range of values: ATsoat = 18 - 180 deg Fy D =) = 22 mm = .158 = 172 in. 

p = 1,0 - 210 atmospheres = 2060 ~ 3090 psia. 

This correlation was plotted for G = 2.ix 196 lom/hr-ft2 for 2000 
psia and 2500 psia and compared with the data of Doroshehuck and Lantsman 
(9) for tubes of bh mm (.158 in.) and 8 mm (.316 in.). At 2000 psia the 
correlation does an almost remarkéble job of predicting CHF for D = 8mm, 
At 2500 psia the correletion is very close to the experimental 6 mn date 
although approximately 20% low at low subcoclings. The correlation 
does not fit the mm (.158 in.) data as well, with differences of 15% - 
30% at both pressures. Figure 9 shows the correlation at 2000 psia for 
the mm diameter. 
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Gambill: (fref: 30) 


i 
i = ¥ = e ee o 
(a/\)er = K hfe pv ("te ali 2) 1 ‘y e Cy Alsat 
v : 








Afe 


+ K? (k/D) Re™ Pr? (t, = ty)go 
where K = 12 ~ ol7, Kt, my & n Gepend on the non-boiling turbulent 
flow cerrelation used. 
Range of values: V = .05 - 17h ft/sec, p = 4.2 = 3900 psia, 
OT.a4 * 0 - 506 deg Fy a = 1 = 57,000 g's 
This correlation is of the superposition form where a subcooled 
peol beiling burnout term is added to a non-boiling forced convection 
term, typically the well know McAcams correlation for turbulent flow. 
Gambill leaves some latitude for the evaluation of the constant 
coeffiecients and exponents. Evaluation of the wall temerature (t,) 
in the forced convection term is doe by an auxiliary plot of AT..4 
V8. Tent/Terit based on Pernath's correlation of wall temperature at 
burnout (16). This is explained in detail in the section on Bernath's 
correlation, The correlation was clotted under different conditions 
with the following results: 
1) D*= .2h2 in, & = 3.69 x 10° lon/ft“hr, p = 90 psia (figure &) 
The correlation fell about 35% lower than the experimental data. 
2) D= .079 ine, G = 3.69 x 10° tbm/hr-£t2, 7.38 x 10° lbn/br-ft*, 
& 1.76 x 10° lon/hr-ft?, p = 90 psia (figure lh): The correlation 
fell as much as 100% lower than the experimental data. 
3) D = ,180 ins, p = 2000 psia, G = 2.1 x 106 lbm/hr-ft* (figure 9) 
The correlation is conservative (low), but does not show the trends 


of the data completely, and is in error by as much as 25%, 
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BPernath: (ref: 16) 
(q/A)er = Rye (ty - by aso 
where: hh, 2 5710 + L8 F p72) Btu/hreft"=deg | 
and: ty),, = 32+ (9/8)/ 57 In(p) = Skip) - y | dag | 
p+ i 


Fance of values: Round Sibes, p = 500 = 3000 osia, VY * hed = 

5); ft/sec, where D is measured in feet, and p in psia, 

This correlation is based on a limiting heat transfer coefficient 
and a correlation of wall temeratures at burnout. The correlation 
was plotted (figure 9) for G = 2.1 x 10° 1bm/nr-ft2 at 2000 pete, and 
D = ,130 in. for comparison with experimental data. The correlation is 
slightly conservative predicting this data within approximately 52. 
Pettis: (ref: 17) 

(q/A)... * 10° x 0,28( h,,/'10 7223 (4 + G/107 )2 et 22L/D  seu/nratt? 

Range of values: G * 0.2 x 10° ~ 8.0 x 106 Lom/hr-£t", p = 1650 - 

2159 psia, Altgat = © - 180 Btu/lbm, L/D = 21 ~ 365. 

This correlation is compared with experimental data on figure 9, 
It ean be seen that it follows the ceeneral trend of experimental data, 
With variations of 15% to 40%. 
Griffith: (ref: 13, 17) 

(a/Mer = £( Pie ) (F) 5 (hg = hy) |! (Py =z) 2" 7 

1 ¢€ 


(F) = 1+ (Fey x 107 ) ~ ou) Ky + 055 x 1077 (Re, Ky 2 


Mp "Yo D fa 3 Ky = 7162 (Toe - TH) 
Mil e “ig : 


and, f( p/p,) is presented graphically in the reference. 
Range of values: p = li.7 - 3000 psia, ¥ = 0 ~ 110 ft/sec, 
LD Tage = 0 = 28 deg Fg 
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The correlation was plotted and compared with sata av two pressures: 
90 psia (figure &) and 2000 psie (figure 9). At 90 psia the correl tio: 
4s outstanding and does an excellent job of predicting the date. 4% 

2000 psia the correlation is higher then the data, yet only by $7. 
“irghak, et al: (ref: 13) 

(q/A)er = 460,000 ( 1 + .0131 p )( 1 * .0050¢ Ae, )(1 + 29365 V) 

Range of values: p = 25 - & psia, ¥=5 -1h5 ft/sec, D = .21 ~ 

eliS ing, where Avent is measured in deg F, and (a/Aor in Sth f/nreft*. 

This ecrrelation was compared with experimental data for .2h2 in. 
diameter tubes at 99 psia (figure ©) and .305 in. diameter tubes at 30 
psia (26). In both cases the correlation curve was conservative (lower 
than actual data), 

DuPont: (ref: 8) 

(@/M op = 595,000 (1 + .OhV)(1 + .0055 Bigab)(1 + 2005S AT sap) 

Fange of values: V = 8.5 - 0 ft/sec, AT... = 40 - 120 der F, 

D= .3 = .05 ins, p = 60 - 1000 psia, Al sep * Subcodling Gradient * 

(Texit - Tintet) per 2 foot length of tube. 

DuPont's correlation is similar to the one developed earlier by 
Mirshak, et al (13); a combination of linear factors of velocity and 
subeooling, in the DuPat tests a two foot length of tube was used, 
Subecoling Gradient (A Tseq)s the temperature difference through the 
tube, was the last factor used to correlate experimental data, using 
the form of a previous correlation for anrmlar geometry. 

Doreshchuck and lantsran:(ref: 9) 

(a/\ ey * (@/A)g # ACD) | (.316/0)+33 i| Btu/hrft* 

where: (q/A),. is determined from a plot of (q/A)o, vs. Ah,., fer 

8 mm (.316 in.) tubes for the pressure in question. 

Range of values: D © .117 ine = 316 in., p = 50 = 170 atm = 7386 - 

2500 psia. A(p) is a tabulated function of pressure. 
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po atm psia A(p) = Ptufhr-ft° 


50 738 2.56 x 19° 
60 1140 2,12 x 10° 
100 1470 2.31 x 19" 
U0, — 2060 2.13 x 10° 
170 ©2509 1.99 x 10° 


Rased on the 8 me (.316 in.) data of Doroshehuck and Lantanen, the 
correlation predicts critical heat flux for ths other diameters with a 
maximum deviation of 30%, for fixed values of the cther variables. “Ath 
other exnerimental data the variation is slightly creater. 

MIT Small Diameter: | 
(q/A)_., = 1.738 x 10° (eS + £01 Ahgat)(6/59)*?? (prey fga) 
where: G, = 10° lbn/br-rt? 
D, @ 0.01 in, 

Range of variables: D = .05 - .305 ins, p = 30 - 90 pela, G= 2x 10° - 

15 x 10° 1mm/nr-ft?, L/D = 5 = 60, Abia, = 0 = 150 Btu/ltm 

This equation was developed by the authors for the low pressure 
data from MIT and fits ‘she majority of the data to 2 25%, 
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RESULTS Or LP. KIMPATAL STUDY 
Scope of Experimental Study 
As mentioned, there is a great deal of expcrimental data available 
in the literature. Although each investigator has attempted tc be 
somewhat systematic in his own methcd of attack, there are very few 
cases where data’ from cne investication can be ecvmpared with that from 
another. Generally there are substantial differences in the variables, 
and lecking a good correlating equaticn, corparison, except in a very 
peneral way, is irpossible. The experimental phase cf tris study was 
conceived with this situation in wind. Ornatski and Kichicin (h, 5 
and 6) have ccnducted tests with tubes cf two millimeter diameter, and 
four mass flow rates, over a range of pressures from 147 psia tc 28h0 
psia. Part of the present study was to complete the low pressure end of 
the 2 mm map. Also, using the same size tube and flow rate, other 
parameters were varied. The results of this investigation, when combined 
with the data from Crnatski and Kichigin, present a fairly complete 
picture of CHF for tubes of two millimeter diameter. 
Experiments were conducted to yield the following sets of datas: 
a. CHF at varicus subcooling for three mass flow 

rates (G = 3.69 x 10°, 7237 & 10°, 1h.76 x 10° 

lbm/hr-ft* - identical to those used by Ornatski 

and Kichigin) with D = 2mm (.079 in.) for exit 

pressures of 30 psia and 90 psia. 

bo. CHF at various tube diameters, with diareters 
varying from 0.0k7 in. to 0.22 in., with mass 
flow rate held constant at 3.59 x 10° lom/hr-ft*, 


and at constant exit subcooling of 35 Btu/lbm 


2h 





Ce 


i. With L/D held constant al 2 vilue 
CP L0, 
ii. With L held ccnstant at a value 
of 3.15 inches 

Ci for various L/D, for the sam: mass flow 
rate and exit subcocling as in part >, above 
with D = 0,079 in. L/D from 3.5 to 50. 
CHF at various subccolings for the sare 
mass flow rate, with L/D of 10 and 
D= 0.2h2 in. This data was obtained in 
order to present a ccmmariscn vith the 
corresponding D = 0.079 data, and es 
substantiation of the CHF verses diareter 
data. 
CHF, for the same flow rate, exit subcooling 
and L/D, 5 = 0,09) in, with changing entrance 
conditions. tn addition to the calming 
length form us: d in the majority of the 
study, four entrance forms were investigated. 
In all cases, the entrance was placed 
imrediately up-stréam of tre heated length. 
Entrance forms investirsatec were: 
1. Sharp Edge Entrance. 


ii, lS-degree Chamfer, 


iii, Rounded Entrance, radius of 


entrance equal to 2 diareters, 
iv. Reentrant, length of reentrant 


section equal to 10 diameters, 
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f. “exifiam Heat Tiix. An inveshies* ich 
of the svster capabilities, using 
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. Pie 


tide 


na. 2 OF tahtaee with 1/9 of 
7.75 and the me imun availablo flow 
rate an? exit subrooling, tc determine 
the raximam obtaimble CW at this 
diameter. 
Results 
A tubal of 117 CHE data pints were ottaines in the exoerimenta!l 
phase of this study. Of these, twelve were obtained usinr the Berylliwn 
Copper tueing. A deteil<cd deseription ci the nroblems enccunters¢ 
waeh thie tubing te corbained 14 Avpendix Jan 
Figures 3 and show the CH data obtained with [ = 0.079 in., 
for the three mecs flow rates. These firures show the sare general 
Gependence of CH on mase flow anid subcoolinge os the majcrity of the 
other experimental data. This lew pressure CHF data “oes show the 
“inimum CHF in the subecolecd region more clearly than does the date at 
higher pressures. [It can be seen thas a minimum CHF is enecuntered at 
approximately 30 Btu/lbm subrocling, and that the subeccling for 
vinimum CHF increases for increasing mass flow rate, CHF beccmes 
nearly independant of subcociing at the lowest flow rate, at 39 psia. 
The dependence of CHF won tube diameter ic shown in figure 5. 
Investigations were ccnducted «ith L/D held constant and with L held 
constant. The crossing cf the two curves merely Shows the dlanoter at 
which the selected values of L and L/D ccincide. The data cbhtained 
with L held censtant shows a greater CHF fcr the larver lianeters, 
reaching a meximum increment of ap roximately 15% at the 0.242 in. 
diameter tube. The value ef CHP can be seen to beccme less depesmient 
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upen diameter fer the larger tubes, Limits of the apparatus precluded 
investigaticn of diameters larger than 0,02 in. at the flow rate 
Selectec far this study, 

The comparison between CH data at D = 9,079 in. ani tnat at 
D = 0,212 in. can be seen in figure % This cata vrcevides an evaluation 
ef the validity ef the curve of CH verses diameter (figure 5) which 
was determined for a particular value of exit subccoling. It can be 
seen that propcrticnal difference in CH between thse two siameicrs 
is nearly constant over the range cf exit subcoclings considered, 

Investicaticn cf the effect ci changing lenrthediameter ratio 
upon CHF (figure 7) shows that CHF increases for decreasing L/D, 
This data shows ccnsiderable scatter at the low values cf L/!. 
The test section construction precedure, discussec in detail in 
Appendix ©, was such that very small fillets cs scldcr remained a 
each em ot tne heated? lenrth, In these small sec*ions, this limited 
the accuracy of determining heated length, It can be seen that Cr 
approeches a ~inimum value as L/P is increased, and that CHF is 
essentially indevendent of L/P tor values of L/D greater than 0, 
This is in agreement with the results of Gergles (1, 11). It should 
be noted that the value of CHF for L/D of 10 (the generally repcerted 
minimum length in the Russian papers) is scme 15% greater than the 
minimum value, 

Results of the investigation of Entrance ifiect are not plotied, 
Folding all flew variablesecnstant and chan,ing milly the fcrm oi the 


entrance tc the heated length yielded the following data: 
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Entrance Form Yeat Flwe 





Sharp iige 3.65 x 10° teu/ft*hr 
LS° Chamfer 3,0 x 10° Rtu/ftohr 
Rounded 3.62 x 10° Bbu/ft2nr 
Reentrant 3.55 & 19° Stu/ft “hr 
%0-D Calming Length 3,82 x 19° Bbu/ftZhr 


The value for the shrrp edge is a mean of Uwe test values, the others are 
Single tests, The 206D Calming Length was whe standard form used for 
the remainder of the exverimental study, the value listed is the result 
of sevérel tests at the sare conditions. Compleie date ior these tests 
are lis’ed under runs no, 33 throurch no, 37 in the ligtinre cf data, 
Appendix C, 

In the stucy cf maximum heet flux, subecclings on the order of 
200 Btu/lb,, and flow rates on the order of 5 x 10° lon/ft*hr were 
ebtained with a test section pressure cf 85 psice. These conditions lead 
toa CH of alrest 6 x 10° Etu/ftehr. Fxaet deta on these tests are 
contained in the listing ci data, Appendix C, runs through numbers 


29%, 30 and 31. 
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IV COMLUSTONS ANT) RECOMNNTATTONS 
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Riscussion of Experimental Fesults 





The data obtained in the experimental ohesw o; this stuay s hows 
the same general denendance of CRF upcn the cther variables as has 
been reported in similar studies. This data is perhaps a little more 
well ordered than that ncrmally reported. This, of ccurse, is 
heartening to the investigators, and is believed tc be the result of 
careful experimentel precedures and constant cognizance of the flow 
Stability rroblem. It is felt that instrument and reacing accuracies 
would ccmbine tc insure an experiment?1] accuracy of five per-cent or 
better. 

The curves of CHF verses exit subcocling terminate at the maximum 
subcooling obtained. This maximum is dependent upcn the test section 
inlet terperature and pressure, The curves show the maximum subcocling 
available consistent with the capscity cf the installed heat exchanrcer 
and the terperature of the city ccoling water (minimum inlet terperature 
= ):7° F), The additional 15 or sc Btu/lbm of subeooling which could have 
been obtained by the installation of a complex test section pre-cooler 
cid net seem tc justify the additional ccmplication, 

The curves of CHF verses exit subcooling show a minimum value of 
CHF, signifying a change of flow regime. This rinimum occurs in the 
subcooled region mainly because of the manner of determining cenditions 
at the pcint of CHF. Fluid ccnditior were determined using standard 
thermodynamic var ables, which are valia for equalibrium states, but 
which de net indicate the physical conditions with inherently none 
equilibrium subcocoled, forced flow boiling, What is needed is a better 
means of describing the physical state, one which would include some 


allowance for the nen-equilibrium conditions as they actually exist, 
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The use of a voicefraction varidlule hat bovr Bugrerves for her ourpa 
(35), but this is very difficalt to avelagte ant oe tepamtent w60 mech 
the same variables as (1%, The use of subtcoliug, as deuercina bv 
a difference in enthalpy cor temereture, ia cconrerigrt, but wusl be 
accespanied by a realizaticn cf its inherent inanility t) dAsecribe 
the actual chys sical phenomenz, 

At the very hich values of exit subroclinc, tir curve mo fo 
seems to become less steep. That is tc ssy that the rate cl change of 
CHF with subcoolinge decreases with incressing subecoling. This iv « very 
slight effect, not at all pronounced at low syster pressures. The 
high pressure date cf TPoreshchuck and lentsman (19) and Crnmtski and 
Kichigin (5, 6) shows this effert more clearly. A possible explenation 
for this lies in analysis of the effect cf subecoling uson the cverall 
boiling mechanism. it appears thet onee high values of subcocling 
are reached, and a large recondensation potential is established, 
further ircreases in subcooling have less efiert 

The inecresse of CFF with decreasing tube diereter is evident in 
figure 5, it has been sugrested that void-frartion inerscases with 
decreasing diameter, This is due tc the manner in which void-frectiom 
is defined, rather than scme vhysical effect, The increase in CHF 
is a direet result of the req:zirement fcr continuitv cf mass flow. The 
size of a vapor bubble is deperdent mainly upon pressure. In sll tubes, 
the bubble occupies a proportionately larger part of the cross section, 
requiring that the Liquid phase assume a higher velecity to maintein 
continuity, It has been repeatedly shown, and can be easily explained, 
that CRF increases with increasing velccity; thus the increase in CHF 
with decreesing diameter. 


The comparison of CH with lencth-diameter retio (figure 7), 
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shows a leveline cif ah a value ef L/* of Aprroxiverte!y ’D, M4? ta 
a resscnabdle mithimum valu thet can me ampblere! 2) experimental “ke, 
where one is atte-pting tt insvre that the elitcte G wie Gi liernnt 
variables ave separated, 

An interesting effect of chanring L/1 can be seer in Lhe plih ef 
CHF verses ciameter (figure 5) The curve for ccnslant tubs leurth 
re lects the effect of simoltanecusly cranring ~iawet=er and Lun whe 
diameter ratic, At the larrer clameter, the Jabe far -cnsbent tube 
length reflects a substantial decrease in L/™, uhus the precter CHF, 
Unfortunately, the limits cf the apraratus were such tha. larger 
diameters cculd not be investigated et the mass flew rate selected for 
this study, It is slitceoether possible that the curve fcr constant 
length weuld show an increase in CH for larger diameters; i.e. for even 
greater reductions in L/D, 

Comparison of CHF date fcr two tube diameters (fizsure 6), shews 
that diameter is one of the rcre imocrtant variables. ‘“Athin the 
experimental accuracy of this study, one can ccnelude that the variation 
cf CHF with diameter is nearly independent of subcooling, et least in the 
Low pressure region, This requires investigaticn at higher system 
pressure, where the effect cf oressu-e becc-esS more prenounced, and 
where hirher levels cf subcoclins are available, 

The data from the entrance effact study shows that entrance reometry is 
unimportant in determining CHF, Here, as throughcut this problem, there 
is a high degree of interaction between the variables. The investigation 
ef changing entrance form in this study was done with a large length- 
diameter ratio tc minimige the L/D effect. As a result, variations in 
flow due te the change in entrance form were likely dammed out at the 


point of tube failure. TDecreasinr the heated length would have intreduced 
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the complication of the If” effect, Tt is safe te ec. nplutle that changing 
entrance 'ore hac nc effcet for Lich valies of L/D. The scmestct lower 

CH fe cie r@atrant turmermy reflect Cie insvadllity due tc compressibility 
in the non-flow spaces around ths reentrams servticn, tut is mre likely due 


%~O @xDe imanval scatter. 


tH 


The maximun@reat fluc daca is ci ioiteresi onlv for future investigations 
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usin, tre partic Gar anparatus emploved in this study. This deta has no 
valve in avalysis cf the cverall CM problem because of the Limitations 
imposed urcn pressure end mass flcw rate by the apogratus, Tests were 
eencuctest to determine the maximum CFF that cculd be cdblained with vorking 


size tube diameter (0.21.2 in.) an? this apparatus, 
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deserived in detai.2 in she secticn in Fartcrs tnfluecrriy 


Beat Flux (Chapter IIT). 
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Bell: (ref: 21) 
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Thé use of SHIS ccrielefion 1! PesvPLeced Rt EweOebres THY Ze Raia. 
Phere is no restriction <1 wné Cute tiam@ter, tut frm other reculits, 
it hos beun shown CHD hae a stron: i twerss tependence or TlanetOrs 
therefcre, this correletion shsuls be ccnsi Jered susefn 
diemeters (0,25 “in, or larrer) only, ‘lthcurh the ari uracy is 15 to 
3% in the case illustrated (figure 9), the cerrleiiocn is relatively 


simple and essy to ap«ly, and is venevally ccnservative in its predicticn, 


Jens & Lottes: (ref: 17) 


hry 


a power funeticn ferns, and woul’ indice-e 


O 


This correlaticr is 
gere CHF at zero subecciine unless a reslistic lewer limit was 
established, Since the correlation wes desirnaed te have ne necetive 
deviation, it is cne of the most ronservative considerec, ts mein 
virtus is its ease of application. 

Ryabev and Berzina: (ref; 36) 

The 2 mm correlation is adequave fcr the Jaus at pressures abeve 
Us:7 psia. it is not consistent with exseriment for the Lower précsure 
data (30 cr 99 psia). There is o minimum sudcocling, epproximitely 
18 Btu/lbr below which the correlation fails tue to Lhe metheratical 
form, 

The cerrelaticn for D ® hmr is adesuate fcr larver tubes, 
(approximately Bin.), but for smaller tubes it suffers frem lack of 
consideration of a jjameter e.fect. Ryabcv and Uergina ccensidered the 
diameter effect in the derivaticn of this ccrrelation by changing the 
leading coefficient, but they intimated that the diameter effect is 
negligible for larger diametors, This ccrrelaticn ceculs be imrcved 


by substituting a function of diameter fcr the leading ccnstant,. 





Zenkevich: (refs 1h) 


Mraluéticn ci tiscosity am ti8 cc Aleta ac T1is MK tonmpyersturt 
proiuces a curve scmewhet like ti. expe~imetty. Fle ale: ic réewseenced 


if thie cerpr@elasice te tc bé used. 

The ebsence af & tube Siem@eter faericr limits tre use. ules of this 
comrelation. THe lower limit of 0.15 in. as toc gmall tc exclude @11 
diameter effect. It is recommenitec that the use ci this ccrrelation 
be limited to tubes cf 0.25 in. in HMaretor cr lamer, where the 


dependence of CHF u cn diameter is ouite s’all. 





Zenkevich ® Subbotin: (ref? 1h) 
This cor-elation, like 60 many cthers, nerlects the effert of tube 


diameter. The lower limit cf 0.16 in, is tco smo1l to neglect this 


bs 


effect entirely. Fer the berrer diameter (.316 in.) data; this 


cerrelation is excellent. The lower limitation on subeceling is sound, 
ti 


awe to the mathenstical fori cf the equaticn. he uyner Livwit has 


nO phretcea!’ valitity ant seems to ce evidence cf cauticn. In spite 
cf the diameter ccnsiderauion, the mathematical expresricn is 
eenvenient, and tie results fcr lorger tubes are quite recd. 
Gambiil: (ref; 30) 

the curve -cr wall terpera ure is suspect, and available cata of 
wall superheat at burnout ices net bear ovt the graph presented. There 
is toc much latitude in chcice of cunstants and exponents te be used 
in the correlation. Frresvumably, if the data is reascnably well ordered, 
the fit of the corm laticn cen be immroved br adjusting the values of 
K, K', m, and nw There are no casily arplied puidelinee for the 
adjustment cf these constarts inciudec in the presentation. This 
correlaticns does include a cinreter tern, but it is enly in the ferced 


contertidn part, and then tc only the (m - 1) rower, where m ie cne of 
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Gambill's "variable" ccns ants, generally on the crder cf %3. This is 
one oi the few correlations which includes small tube diameters, yet 
clearly the correlated effect of diameter is less than that observed in 
the data of ergles (1), Poroshchuck and Lantsman (9), end figure h. 
The correlation is no better than that cf “irshak (13) at lew pressures 
and is considerably clumsier to handle. At high nressures (figure 9) 
this correlation is no better cr worse than many cthers. 

Bernath: (ref: 16) 

For a given pressure, tube diameter and flow rate, tris correlation 
is linear with subccoling, This is a faitly good first avproximation 
and typical cf many correlations. This ccorrelation is rather di: ficult 
to apply, requiring auxiliary calculations to determine heat transfer 
coefficients and wall temperatures. Cf the seven high pressure 
correlations plotted on figure 3, this is one of the best, precicting 
CHF within 5% for the selectec conditions. 

Bettis: (ref: 17) 

The use of this correlation is restricted to pressures near 
2000 psia. There is no restriction on a minimum diameter; therefore, 
the correlaticn will yield low predicted values of CHI for diameters 
less than .25 in. 

Griffith: (ref: 13, 17) 

This correlation was ogiginally determined for several fluids 
ineluding water and it is applicable throughout the subcooled boiling 
region and into net steam generation. Generally, over a wide range 
of variables, this correlation is the most ccnsistently accurate one 
available. Its disadvantage lies in its extremely cumberscme and 
unwieldy form. Computation can be simplified screwhat by using the 


curves presented in WAPD-188 (17) for several cf the recurring variable groups. 
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Mirshak, et al: (refs 13) 


Tre for® cf the ccrrelstepr is liwear 1) @i.beicclim® am!) és 
exhinits the approximate ventenc#® ~f th exnerMiente, dete. ithe 


the ranpe c.f dianeters icr with The ccrreWohun ‘s Mitel, #he 
deperdence of Ci cn @i@meter is negliecte’. The ocllesine pepere! 
%, 
ecments are anpropriate: 
) The correlation ives a ccnservative in’icaticon 
1) The correlation ives a ecnservative invication 
of CHF, but errs by 25% - 50% 
2) There is no censideration of the variaticr of 
CHY with diameter 
3.) The correlation is safe, suick on’ easy te azunlv 
TMPont: (refs 8) 
This correlation was des§enet f°~ correlating wage 'cwedata for 


tubes of a especisic length, Hew te ev-Luate Subecoline Gratient for 


cther tuoe lengths is unresolved. ‘Subeccling gradient is not an 


* o 


independent variable in Lhe equation, beins tied to Heat } lux by the 
Yirst Law of Thermcedvnamics, “ercter, ond flor rate, No diareter 
effect is nected between 9.39 - 0.38 in. This is a re-senabdle 
approximation. Variaticn in pressure between 60 -1000 nsia is 
considerable and shculd be taken iritce account, Since the enueticn 
was derived frem a small sample ci data and fer tubes of only one 


length, its se is nct reccmeenied. 


Poroshchuck & Lantsman: (refs 9) 





This is cne ci the few stutvies which hes been devcted exclusively 
to the diameter efiect oreblem, The reneral treni is correct, but 


comearing this ccorrelatine enquaticn with exoerinmental data - even that 


fror the paper in which this cor eletion is nresented, shows en error of 


to 30%. Part’of the difficulty is due tc the int@rarticn of the variables. 
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Cof¥oarinrg tlis cchewe ®ith ener cemmeagemiel ceta, nqrably that 
of Cipmatski.end Zichtein (1, 5, *,) amd SPr@les (1, 11,) Ble shewed 
the correlaticn tc e@eecd a 20M errcr., “8L Mes © on Wire 4 Se cf 
A [(0/n, yi oe x fit the Ornotski and Kichisin data (i, $ and 4) 
quite well, cut this failed in c@iperisc i with gat: free cthéde sOidics, 

Tt 18 nct weccrememded that this correletic4 o5 user tc fred ct 
CHF; hawemer, it will serve tc annly «a cieme~e~ cametigeraticn tc ctlier 


correlations, within rather bread limits. 








Th’ s eguaticn was develone! act 


uy 


C MurWicen a oredicwen of ST but 
$@ nrevice a feirly ~encise ethemetical Gepresricn fcr the lek oressune 
MET Fata. “Yo prescure term is i> lu@edms the variation of CHF with 
pressure is neither well ordere? nor of lerc-e magnitude. The change 

in CHY with -ressure cver this range is mmich less thar the error of 

tha$ Gaueticn due te the assumpticr cf CMF as e Linear function of 

weve ing, CREF is, in fact, not Dinear with sunecclingy this ferm 


a 


has been assumec merely fcr ecnvenience,§ Tt 


tte 


s recormend«d that the 
use of this enuation be Limite? te deter-ination cf *neral levels, 


ang that the artual data be used in anv -ttemots at accurate prediction of 


CH. 
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PROC PURE BCR OGh 
Firures 3, 4, 10, 11 and 12 provide 4 rranhical presentation of 
CHF fcr pressures from 30 psia to 2040 psia cver e range of mase flow 
rates and exit sabccoolings fer tubes of two mm diameter. These curves 
in conjuncticn with the data of CHF dependence upon ture ciamerier 
provide a means of geraphically detervining CHF for tubes un tc onee 
fourth inch diameter, flow rates from 2 x 10° 1b_/ft? hr. to 22 lbp/ft? hr, 
at any pressure and exit subco ling. A precedure of linear and graphical 
interpolation is recommended for predictine CHF for any reasonable 
small diameter application. In chases where linear interpclaticn is 
recommended , auxiliary curves, have shown this technique tc be 
satisfactory. The accuracy of this eraphical prediction technique 
will exceed that of the few correlating equations applicable to stall 
tubes, and in the majority cf the cases tested, was within 15% cf the 
experimental data considered, 
The following procedure is recommended for predicting CH: 
l. Select desired system pressure, macs flow 
rate, tube diameter and Length. 
2. locate the plots for pressures bracketing 
the desired pressure, 
3. At the selected exit subcocling, read 
CH for the three or four mass flow rates 
presented, for each pressure, 
li. Prepare interpclation curves with coordinates 
of CHP and Mass Flow Rate, and determine CH 


corresponding to the desired flow rate at each pressure, 
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CS. Usine the twc values thus ol-tained, Hnear 
interoolation between th* Wressures unile 
censideration will sive a value ef Cm fer 
the cesired pressure, flow rate and 
subecoling, for a diameter of two mm, 

6. Reter to figure 5, end using the curve for 
ecnstant L/D. deterrine the value of CFF 
for the “iameter in questicn. The ratio 
of this last value, tc that at two mm 
(a value cf CEF cf b.15 x 10° Btu/ft] hr. 
at 2 mm diameter will vrovide satisfac cry 
results) multiplied’ by the previously 
determined value cf CHF, will yield a 
predicted value of CHF at the desired 

conditions cf pressure, flow rate, subcocling 


and tube diameter. 


7. In cases where the desired L/D is les, than 
20, the same sort of ratio correction should 
be performed, using figure 7. for values of 
L/D of 20 or greater, the error in CHF will 
be tec small tc warrant consideration, 

For a rough, quick estimate cf CHE, visual interpolation for mss 
flow rate on the plot fcr pressure nearest the desired vressure can be 
easily accomplished. The diameter correcticn should be included in this 
abbreviated procedure. This same sort of apnroach is useful in 


estimating the effect of gross change inthe impcertam variatles,. 
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Figure (19) - “ependenre of Critical Heat Flu: on Subcooling 
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Figure (12) = Pevendence ef Critical Feat Tlux on Subreocling 
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Experimontatdon wes ccnductec’ an exicetA epparatils wahch ad 
been ccnstructed for a previous suaty (1), en? which bac been used 


Cm 


for several other toiline and Yivéstiravion®, Univ mincr 


modi fication wh 6% Ret for this investivation, The apparatus 

is shown schematically in fipum 1. The fiuid loop con:ists of a 
closed rain circuletion lcop, with a by-pass psrailcl vest section 
line. All components are of corrcsicn resistaom materiels, The 
system contains a mein circulatirre pump, ac’ umutlatcr, the test seetion 
line with flowmeter, preheater, test secticn wit: instrumentation, 
main loop by-pass line and heat exchancer with city water cocline. 
Auxiliary to this are a fill purp, sunply tank, cesassing tank anc a 
contimuously operating derineralizer, Power is supplied tw the test 
section by a pair of a.ce motor - d.c. generavor sets, 

The main suoply pum is a two-stage turbine wype pump, driven by 
& 3 hp iniucticn wetor. The purp ic capzbie of Aslivering 3.6 callens 
per minute at 260 psi. A pladder type accumulator located at she 
exit side of the pump serves to “amp out pressure fluctuations. Flug 
valves at each en? of the byepass line previde contrcol for the test 
section pressure. Test section flow rate is determined by a Fisher 
Porter flowneter ard controlled vith a Heke needle valve. This needle 
valve elsc provides the 150 psi up=strearm precsure drcep necesssryv tc 
insure flow stabilitv. 

The metcr-generstor sets consist of Lio v., 3 phace synchronous 
motors anc 36 kilowatt d.c. generators, comected in series, ech 
generatcr is nominally rated at 12 volt and 3000 arps. 4 totai power 
of approximately 60 kilowatt is available at the test section. 
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Instrumentation is avallacle icr menitcring pressure levels, 
temperatures, test section flew rate and test secticn current end ¥cltege. 
Test section exit pressures sre reac cn bourdcn type pretcsure gauges 
with specified accuracy cf 0.5 psi or better. The gauges are calibrated 
pericdically with a dead weight tester, Tcst section flew rate is 
determined with t®e Flowrater “eters; the range of flow rate is varied 
by changing tubes ani floats. Meters are calibrated in place and 
calibration curves are available. Test secticn power is obtained from 
the test section voltege and current. Voltare acrcss the test section 
is read directly with a destcn variable scale, calibrated d.c, voltmeter, 
Current is read across a “ational Bureau of Standards shunt, with a 
specified calibration of 60.17 amps/nv. 

All temperatures are determined by thermcccuples manufactured 
from 30-gauge, duplex, Ccpper-Constantan thermecouple wire. Test 
section inlet and exit temperatures are measured on thermcrouples 
inserted intc the flcw with Conax fittings. Output voltages of the 
thermocouples and current shunt are displayed upon 2 Minneapolise 
Honneywell Brown, O = 26 my., single channel recorder. 

A much more detailed description of the apparatus is contained in 
ref. 13; however, certain modifications have been made, The more important 
of these include: 

ae The braided copper flexible power connection 
anc the rigid upestream jumper from the 
test panel to test section bus connection 
were replaced during this study by Mackworth 
Reese Flexible Dry Jumpers. This requires 
that the test secticn be suvported. Support 
is pervided by a laboratory support and clemp, 
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Ce 


with electrical irsui Ass, 
The braided copper nccer sumply cebles lave 
been replaco? .ith Vackwerth “esse, water 
coolcd cables, These cahies are erouled by 
city water tanped off ci the noe. exchayeer 

%, 
eoclant supply line, discharged bo a drain. 
The preheater has been vod’ fied by the 
addition of three 5 Kilowatt imrersion heaters, 


with ceparate eccntrolis, These heaters are in 


sertées with the orisina 6 tllovatt var’able 


contro] heater and previde any desired 
preheating frem 0 te 21 kilowatt, with 
steady state operation, 


The steam heater in the supply tank hee 


obtained fror the cegas tank heaters alcne. 
The thermeccounle selection switch now 
consists of twc Leeds and Northrup, 12 
positicn selectcr switches, ins ries, 
Thernccouple wires eré connected through 


crew type terminal boards, thus recucins 


ee) 


set-up tire, This rodificaticn was made 
in ofder to incresse the number cf 
thermeccuples which could be natched 

to the rercrdcr, Hesictanee bests showed 


resistance as a 


hfe 
> 


mo aprreciable change 


result of the added switch, In this study, 
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oniy the primery seit? Wis used, 
When using ule ver’ Shert test secticns, 
er theme cit tyvaee™*" ~ ci, che test 
secszicn retistance was so lew ay t. 
wake tie mctoreronergetcor “*oogrcl 

%, 
highly erratic. Te alleviate this 
prcbler, a seriés resistenec was 
fabricated. Thies resisvance consisted 
of a 10 inch length of 9.2h2 inch type 
39): shainless steel tubing, with ous 
ecnnections, place? Slectrieelly in 
series with the Th whtalen ait cociad 
br etty water. This giiat reeictennaé 
eauced an aiditicnel ecltage tro” of 
aprreximetely 5 velts aercss *he norer 
Supply, which irmreved the mec. set 
contre] considerably, Power to the test 
secticn was cetervired by reariinre rottare 
drop across the teét section only, and by 
reading current threurh the series 
combi naticn, 
In the tests where very hich pewer 
applicaticns ware required te reach CH, 
the test section exit bus connecticn was 
fcund to heat uy excessively, causing 
seoking and in ne case 2a small fire, rot 
tc mention creat ir-onventence when chencing 


test sections, In view of this, e bus 


comnecticn coolant scheme was devised, A 
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teinch copper ninple in’ twc el ows wer 
silver brazed tc the brass bus connecticn. 
City water, fror the nower lead cooling 
system, wes diverted throurh this ninple, 
drawing hee. from the bus conrection, Heat 
balance checks when cperating, comparing 
power input to fluic ter-mersture rise thrceugh 
the test section showe’? that tris ccoling 
device hat no e’ fect on the test date, No 
further overheating problems were enccuntered 
once this ccoler was installed, end it was 
used whenever there was a requivevent fcr 
high power input, 

Test Section 

Pecause each test destroyed the test section, ccnstruction was kept 
&s simple as possible. 4 test section wes constructed by facing a 3/h inch 
brass pipe cap and drilling a hole in the center of it. 4 secticn of 
the test tubing was silver soldered intc this hole. At the other end 
of the desired heated length, a 5/16 inch brass bushing, about 5/8 inch 
long, drilled to the 0.D cf the tubing was silver scldered to the 
tubing. Heated lengsth was measured from fillet to fillet. A calming 
length cf approximately 29 diameters extended beyond the up-stream 
bushing. 

Sections were installed in the apparatus bv threadine the pipe cap, 
the down-stream end, to the dcown-stream plenum chamber which contained the 
exit pressure tap and exit thermocouple. The up-stream end was installed 
by inserting the calming leneth into a Conax fitting sealed with a 


neoprene washer. brass bus cennections were installed on the bushing and 
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adequate immundarce Mavcice. 2G atl 
lieting of Bate (Eocenbge , rate GO, 30 we ad). 
%, 

€ ‘ iti | ~i y 
Beryllium-Coppver Tic. 

At the besinirg ef: We cate timentaticn, tain df Of ello, Me-'u 

2 s ee ae ‘ * ms fn \ 

wee orceze], with immite Chomeber ef 9.089 berrec (7° man}. Ttewaeée f@1t 
that this tubing eMol? be @eeier th fem FLY) “han bee i vce 3% stainless 


oi i - : * ae ~ , = 7 . - 
siyee. whieh he. Saw ce: Ln ure come Stadio. gvene! Lest s¥etchens 


~ 2c o - oe s = 
were mei uo ani tecve were sun. “Tetertinctol), coumidenmable difficubty 
Was Oreecinterec using the Ueefu Ame nl Come tle =r) temcmant 


proplems are discussed belo. 
About the only iabricevicn efvantars that dho Se-"u tulling did we 
Wee it ceulr be om} _nto desired ilengis bedi, Untorbunetaly, the 


pt 7 a + yah ~ % & . 
Seaey Was not siatie at teMmtratures in exress of 11° degrems F. 


C 


MM test secticns were constructed by sidyver scléeriag, the tutine near 


Gee Scliered cejiaections wiervcnt 2 phase change amd this was apcommanied 


q 


by a change in electrical resistence ef almest YO, Since this change 


in resistenre occvred only in the repion where heat was applied for 


ae 


anea eoysed an uneven arxia? 


t-te 
{n 


sceidering, the lecal charce in res 
Mwetribwbiem cf heat flux. "late coteained hw maaserine tote). test 
section voltare drop anc current did net incicate the actual SE. 
Test sections were conscructed using seft seldery but the hirh 
emit temperatures efitcurtered durin> overetion cm iBed the eeft sclder 
to cuddle *ust befcre roachine CIE, “ec leaks develcped, Sut thie was 


felt to be less than satisfactcry, Ascher temerature (99619) soft 


solder wac tried, with nc beter results. 


65 


Se = 
~~ «© -— 
a ' ae 
® —-— -— emt) 
~—— elie 
—_ ‘* = 7 
= * - ees: eee — 
— a ' ‘ ~ nm » 4) oe 
a — tA Deas 
ee Fm op, 
= ; — nee doce) 
= «te onde 2 
—< = & - 


ei ——_ § - Ja, 


OE) NO = ee! a «to mr 


we 
>» & 













o tx a. °°. © |) aoa 1 


4 


—e = a 
er ee) ed Pee Lee © ile: ee 


SE ee ee Sake tec Ss 
SF ey ED Geer cee eee — = 
en Lu. © © pes a 
eye. Gk woes cet age 
ite ee ee 
OO Dom: ee ee oe A ilies tee Or encore betes 
eS! eet ee cake te | 


‘~~ 
ue 


7 


7 








>on 













The lower resistance Ee-Cu tubins csused mp. set ccntrcol problers, 
The series resistarce previously mentioned alleviated this sroblem, but 
tended to cerplicate the procedure. 

Cll date was obtained using tis ‘e-Cu tubing. “ubsequent ccrparison 
with similar cata cbtained with tvpe 30) stainless steel tubing shewed 
that the Be=Cu data was inconsistent, for either the soft scldere? or 
Silver soldered sccticns,. 

A corrcsion check wes meade with the Bee-Cu tubing, After one hour 
at a heat flux of atcut one-third of CH, the inside of the tube showed 
distinct blecxeninr, using derassed and cemineralized water, 

In view cf the various difficulties, it was conrluded that the 
Pe-Cu tubing was unsatisfectcry for experiments of this type, and thick 
wall, type 30) stainless steel tubing, with inside diameter cf 0.079 
inches (2 mm) was obtained. 

The stainless steel tubine could be silver scldered satisfactorily 
and the data obtained from it showed considerably rore consistency. 

The data presented as the result cf this investic¢ation was obtained using 
type 30 stainless steel, except where noted te be Be-fu (runs £2 through 


#13), and except for the D = 0,180 inch tubing which wes type "A" Mickel. 
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th@ apparatus, 220 mir wer cléa Trae dhe s-stew, dodth under the Matic 
head of thé sevgs wenh, ~n@ wice tas bam oir aleting water threurh 
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fee 
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panee 1b Wa ceslirer to Hrevemnt @Ble al CA ag & parbiculer Gxt 


susccoliag, it was necescery ic d@vcrine anc estaclish a cesirec Cest 
Section exit Lemperstire. This yee cerplicated wy Uhe .act. that eit 
Lemmereture is fependent uten tlc tris) Ltmoerealiire, which cculd ve 
controlled by the preneater, an? Uron tile nest i.iwx in the test section 


In order tc obtain the desired exit te-pcrature and »re#sure at th 


inc pease the “Niowar to the test sactien 


Sa 
t- 


rs 


Gint of “iT, it war necessery 
Very slowly, an’ t cconbinaally adjust the ted. seeticn preisure &nd 
the preheater controll ny the inlet terperature. In sou eave, the 

el Kilowéttsof preheat eveilshle war nov suflicient, in whith case the 


flow of ccoling wa:er to the hea. excraltiyer was reducel, This. preeedure 


was difficult in that wre exit telmers ure showed a long (about 2° minutes) 


CG 


response tire to aijausiments in the flow &. mity water cculing time flew 


threush the heal exchancer. 
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Often the mrocess of obtainin tne tesired exit subecoling es 6 
very delicate one, This was nperticulary true when ‘ie cure of Clb 
versus subcooling was approximately parallel to the curve cr heat flux 
versus subcocling as ‘etermined by the rirst Law of Thermeudynamics, 

In the large majority of the tes= runs, it was possible ti obtaiz 
the desired valus. of exit subcooling within 5% at the noint of CHF. 

In these few cases where exit subcooling varied mcre than 5% frcm the 

desired value, an auxiliary plot ci CHF versus subcocling was mace, to 
determine the acpropriate value, or to determine that the variation of 
CHF with subcooline wes not sisnificant in the range cf subcooling in 

question. 

For the purpcse of this study, CHF was defined as the heat flux at 
tube failure, While the heat flux thus determined my exceed that at 
the departure from nucleate beiling (DNB), heat flux sc defined should 
coincide closely with the heat flux at the pvoint of transition from 
nucleate to film bciling. Wo burnout protecticn device was used in 
this study, as it wss felt that less time and materials were involved 
in the fabrication of a new test section fcr each run than in the 
construction cf such a device. Furthermore, there crn te no dcubt 
that a genuine meximum heat flux has been established when the tube 
ruptures, 

All experiments were corducted with the test section in the 
horizontal pcsition. System water was degassed only periodically. 

Care was taken to insure that no air was trapped in the system when 
replacing sections and extensive air bleeding was done before commencing 
anew run, The inside surface of the test section was given a thorough 
cleaning with acetone pricr to installaticn. No additives were 


introduced into the water, 
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In almost all cases, UVhe tue ruptured within three ciymevers 
of the exit ccmection. Thus conditions in the exit pletu chamber secre 
considered to de representative of those at she ocint <i failure. A 
recent study by Lopina (3) hes shown that the pressure within the vilenum 
is virtually identical with that at the downstream end cf the hanted 
section throughout the subccoled region, The test secticn and plenun 
were wrapped with fiberclass insuleticn to greater than the critical radius 
to minimize heat loss, 

Because the recorder used presents cnly one channel, it was 
necessary to estimate certain quantities at the point of CHF; however, 
near failure, very small increments of vower were anplied, and all 
quantities were checked frequently. Lata not obtained right at failure 
eculd be readily inferred from that recorded just prior to reaching 
CHF, All changes in conditions were made very gradually during « run, 
and ence the procedure was well established, it was possible to obtain 


a set cf data, e.g. one cata point, in eanproximately forty-five minutes. 
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APP™NDIX F - SAMPLE DATA SHEET 


OCS adn. 2 65 psia Run # 4S 
S15 in. tf 40 | pate 2/25 
S.69 x 406 Lom/pe2 py a w= 126.5 Vo0/py 
7 74 % on Flowrater FP-2-21-G- lO/ac 


QO rr Yo 


Flow @ Tin mv , Toyt mv, p — _ Ee mv Variac % 
C a4 °r> | C \%9 °F) _ 
O02 | 2.70 anc 
- 13.04 750|9.06| -— | 
v - oe E B00 9.66) & 
= ; 
- la20 | 76.19.08 eee a 


ef = aaa [78 [936 1090 [ 


At Failure 
[= [ose [sea | 75 oss izes] — 


S2°r 2GO0°F &9.7 psia 





Heat Balance 


ALT= (25,.5«92 = 12500 


= 


B°1(60.17)(3.41) = 6.O5*«B52*6017*34l = (Z1ISO 
OK 
hsat £90.85 
chr = £'1(60.17)(3.41) h bulk 228.6 
TD *L/y 44 7 2 
Neat GC! BM/1om 
Cur = 10.5512 .6S*( 60.17) ( 541 ) - 4.22 Btu/ p42 an 


THO .O19316 //5 dd 
Obtain Temperature from Thermocouple Conversion Chart 


nh sat = hatp h bulk = h at Toyt 
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APPENDIX G 


DEVELOPMENT OF THE MIT SMALL DIAMETER CORRELATION 


This equation was developed by the authors for the low pressure 
MIT data and fits the majority of it to = 25%, 

(Q/A)op * 14738 x 10° (L.L8 + .OLAKgat) (G/GQ)*9 (0/04) (fa) 

The correlation was developed in the following manner: The 30 
and 90 psia plots (figures 3 &% ) for 2 mm data (.079 in.) were linear- 
ized to determine the slope and the intercept of the part of the equation 
dependent on subcoolong. The flow rate dependence was determined by 
plotting on log-log paper, cross=plots of the 2 mm data at constant 
subcoolings. Also plotted were cross-plots of Bergles'(11) and 
Wessel's (26) data for constant subcoolings. A mean value of the slopes 
was used to determine the exponent of the flow rate term. The same 
technique was used for diameter and L/D dependence, using the data of 
figures 5 and 7, and the data of Bergles (11). The leading constant 
was determined by the zero subcooling intercepts and the values of 
the constants G, and Doe ‘lo pressure term is included as the variation 
of CHF with pressure over this range is neither well ordered nor of 
large magnitude, 

The correlation was plotted over the MIT data and small adjustments 
were made in the constants and exponents to get a better all around fit 
and balance out the positive and negative deviations. 

A further refinement of the correlation would be to make the 
heat flux a weak 2m or 3rd order dependence on subcooling, eliminating 
the linear approximation, A minimum could be introduced at about 
Ahsat = 20 to 30 Btu/lbm, 
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APPENDIX H_ = NOMENCLATURE 


a - Acceleration in g's 

A ~ Heat Transfer Area, £t¢ 

Cc = Specific Heat, Btu/lomeder F 
D - Tube Diapeter, inches 


«- Acceleration of sravity, £t/see" 
~ 32.2 lbpeft/lbm=see® 
~ Heat Transfer Coefficient, Btu/hr-ft°~deg F 


Enthalpy, Btu/lbn 


Mechanical Equivalent of Heat, 778 ftelbs/Btu 


Thermal Conductivity, Btu/hr-ftedeg F 


Pressure, psia 


es /- > |S = 
I 


Prandtl Number 


qg/A « Heat Flux, Btu/hr-ft¢ 

t - Temperature, deg F 

? - Temperature, degrees Fankine 
ATaat - (Tsat - Tb)» deg F, Subcooling 


Vv ~ Specific Volume, ft? /Lbm 

v - Velecity, ft/sec 

v ~ Viscosity, £t°/hr 

-- ~ Viscosity, lbom/hr-ft 

ep ~ Density, lbm/ft- 

c - Surface Tension, 1b ,/ft 
ga 


Ah ° (figng - hy), Btu/lbm, Subcooling 


"2 





Subseripts 


bo # ft the point of burnout or critical heat flux 
er - At the point of burnout or eyvitical heat flux 
b @ Bulk or mixed mean comiitiows 
crit. - Critical” 
fg ~ Saturated Conditica e Varor Condition 
» Vapor Properties 
1 » liquid Preperties 
pb - Fool Boiling 
8 » Saturated Conditions 
sat - Saturaed Conditions 
Vv - Vapor Properties 


w « Yall 
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wy 
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